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ABSTRACT 

We perform hydrodynamical calculations of core-collapse supernovae (SNe) with low explosion en- 
ergies. These SNe do not have enough energy to eject the whole progenitor and most of the progenitor 
falls back to the central remnant. We show that such fallback SNe can have a variety of light curves 
(LCs) but their photospheric velocities can only have some limited values with lower limits. We also 
perform calculations of nucleosynthesis and LCs of several fallback SN model, and find that a fallback 
SN from the progenitor with a main-sequence mass of 13 Mq can account for the properties of the 
peculiar Type la supernova SN 2008ha. The kinetic energy and ejecta mass of the model are 1.2 x 10''® 
erg and 0.074 Mq, respectively, and the ejected ^^Ni mass is 0.003 M©. Thus, SN 2008ha can be a 
core-collapse SN with a large amount of fallback. We also suggest that SN 2008ha could have been 
accompanied with long gamma-ray bursts and long gamma-ray bursts without associated SNe may be 
accompanied with very faint SNe with significant amount of fallback which are similar to SN 2008ha. 
Subject headings: gamma-ray burst: general - supernovae: general - supernovae: individual (SN 
2008ha) 



1. INTRODUCTION 

A massive star with main-sequence mass above ~ 10 
Mq is thought to end its life as a supernova (SN) after 
forming an Fe core at its center. The SN is triggered by 
the gravitational collapse of the Fe core, thus being called 
a core-collapse SN. The mechanism that leads to the final 
emergence of an SN from the collapse is still under debate 
but observations show that ejecta of the SN normally 
has a kinetic energy of ~ 10*^^ erg. Currently, however, 
theoretical attempts to simulate the whole explosion of a 
core-collapse SN have not obtained explosion energy as 
large as 10^^ erg (see, e.g., Janka et al. 2007; Bruenn et 
al. 2009; Burrows et al. 2007; Suwa et al. 2009). 

If the explosion energy is low, the inner part of the 
star falls back onto the central remnant and only the 
outer part of the star overcomes the gravitational po- 
tential. The idea of the fallback was first introduced by 
Colgate (1971) and many studies have since investigated 
the effects of the fallback onto the central remnant, e.g., 
black hole formation [e.g., Chevalier 1989; Woosley & 
Weaver 1995; Fryer 1999; MacFadyen et al. 2001; Zhang 
et al. 2008). Recently, more attention has been paid to 
the outer part, which eventually escapes from fallback, 
thus being ejected. The ejecta might be observed as an 
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SN {e.g., Fryer et al. 2007, 2009) and could produce 
the peculiar chemical abundance patterns of extremely 
metal-poor stars {e.g., Iwamoto et al. 2005). As this 
ejecta has a kinetic energy just above the value required 
to overcome the gravitational potential, it is expected 
to have very low energy. If enough amount of ^^Ni is 
also ejected or the ejecta interacts with the circumstellar 
medium (Fryer et al. 2009), this ejecta might be observed 
as an SN having very low line velocities. 

In this connection, the peculiar SN 2008ha is a suit- 
able object, with sufficient observational data that can be 
compared with the fallback SN models. SN 2008ha was 
discovered on 2008 November UT 7.17 (Puckett, Moore, 
Newton, & Orff 2008) and was found to be one of the 
faintest SNe ever discovered (Valenti et al. 2009, here- 
after V09; Foley et al. 2009, hereafter F09). It was found 
in an irregular galaxy UGC 12682 at a distance modulus 
of /i = 31.64 mag (F09). Adopting a galactic extinc- 
tion of E{B — V) = 0.08 mag and little host extinction, 
the peak absolute F-band magnitude was found to be 
as faint as —14.21 ± 0.15 mag (F09). From its spectral 
similarities to SN 2002cx-like Type la SNe, SN 2008ha 
was classified as a peculiar Type la SN. 

SN 2002cx-like SNe make a class of peculiar Type la 
SNe (see, e.g., the supplementary information of V09; 
SN 2002cx (Li et al. 2003; Jha et al. 2006) and SN 
2005hk (Phillips et al. 2007; Sahu et al. 2008) are well- 
studied examples of this class). Their spectra do not 
have strong absorptions of Si and S at early epochs ^, 
which are the characteristic features of normal Type la 
SNe. Line velocities of SN 2002cx-like SNe are very low 
compared with normal Type la SNe (Branch et al. 2004) . 
They also have peculiar light curves (LCs), which decline 
slowly in spite of their low maximum luminosities and do 

^ However, the earliest observed spectrum of SN 2008ha before 
the maximum luminosity showed these features (Foley et al. 2010), 
although these features disappeared soon {V09; F09). 
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not show a second peak which appears in the / and R 
band LCs of normal Type la SNe. 

SN 2008ha has additional peculiarities. The rise time 
of SN 2008ha is faster than that of normal Type la SNe 
and the decline of the LC after the maximum is very 
rapid: Ami5(B) = 2.17 ±0.02 mag (F09)^°. Line veloc- 
ities of SN 2008ha are as low as ~ 2, 000 km around 
the maximum brightness (V09; F09). Thus, the ejecta 
is expected to have very low energy. V09 suggested that 
the ejecta mass is Mej = 0.1 — 0.5 Mq and the kinetic 
energy is Ei,in = (1 - 5) X 10^9 erg, while F09 estimated 
Mcj = 0.15 Mq and E^^i^ = 2.3 x lO''* erg. The estimated 
mass of the ejected ^^Ni is also as small as (3 — 5) x 10"'^ 
Mq (V09) and (3.0 ± 0.9) x 10"^ Mq (F09). 

Given the low energy and the small mass of the ejecta 
as estimated from its spectral features and LC shape as 
well as its star-forming host galaxy, V09 concluded that 
SN 2008ha is not a thermonuclear explosion but a core- 
collapse SN with fallback. F09 also pointed out the pos- 
sibility of the core-collapse origin but did not exclude the 
possibility of a thermonuclear explosion. Indeed, based 
on the earliest spectrum observed, Foley et al. (2010) 
suggested that SN 2008ha is related to a thermonuclear 
explosion. Alternatively, Pumo et al. (2009) related SN 
2008ha to an electron capture SN (Nomoto 1984). 

In this paper, we show that the properties of SN 2008ha 
can be explained well by a fallback SN model. We 
first perform numerical calculations of hydrodynamics 
and nucleosynthesis for several progenitor models. Then, 
we perform radiative transfer calculations to obtain the 
bolometric LCs and photospheric velocities of these mod- 
els to compare them with the observations of SN 2008ha. 

In Section 2, we introduce the pre-SN models. Meth- 
ods used in our calculations of hydrodynamics, nucle- 
osynthesis, and bolometric LCs are described in Section 
3. We show our results of hydrodynamical calculations 
in Section 4. The results are compared with the observed 
bolometric LC and photospheric velocity of SN 2008ha 
in Section 5. Discussion and conclusions are given in 
Section 6. 

2. PRESUPERNOVA MODELS 

The spectra of SN 2008ha do not show hydrogen lines, 
which implies the progenitor of SN 2008ha has lost its 
H-rich envelope before the explosion. In addition, the 
identification of He lines in the spectra of SN 2008ha 
is very difiicult because the lines are overcrowded (see 
the supplementary Figure 4 of V09 and Figure 8 of F09 
for line identifications). Thus, both a helium star and a 
carbon + oxygen (CO) star are possible candidates for 
the progenitor of SN 2008ha. 

We use pre-SN models of solar metallicity with main- 
sequence masses of 13 Mq, 25 Mq, and 40 M© calculated 
by Umeda & Nomoto (2002, 2005). As these models have 
a H-rich envelope, a He star is constructed by assuming 
that the whole H-rich envelope is lost either by stellar 
wind or Roche-lobe overflow of a close binary, and only 
the He core remains at the pre-SN stage. The boundary 
between the He core and the H-rich envelope is assumed 
to be at the location X(H)= 0.1 (hereafter, A(M) de- 
notes the mass fraction of the element M). We adopt the 

Ami5(_B) is the decline of the B band magnitude in 15 days 
since the S-band maximum. 
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Fig. 1. — Density structure of the progenitor models. The line 
of p oc r"-^ is also shown for comparison. 

25 Mq and 40 Mq models to construct the He star mod- 
els 25He and 40He, respectively. The He core masses of 
25He and 40He are 7.0 Mq and 15 Mq, respectively. 

The CO star models are constructed by assuming that 
both the H-rich and He envelopes are lost and a CO core 
remains. The boundary between the He envelope and 
the CO core is set at X(Re)— 0.1. We construct CO 
star models, 13C0, 25CO, and 40CO, from the 13 Mq, 
25 Mq, and 40 Mq models, respectively. The CO core 
masses of 13C0, 25CO and 40CO are 2.7 Mq, 5.7 Mq, 
and 14 Mq, respectively. The density structures of all 
the models used in this paper are shown in Figure [1] 

3. METHODS 
3.1. Hydrodynamics and Nucleosynthesis 

Calculations of hydrodynamics and nucleosynthesis are 
performed by using a spherical Lagrangian hydrody- 
namic code with a piecewise parabolic method (Colella 
& Woodward 1984). The calculation of explosive nu- 
cleosynthesis is coupled with hydrodynamics and the 
adopted reaction network includes 13 a-particles, i.e., 
4He, 12c, 1^0, 20Ne, ^''Mg, ^^Si, ^^S, ^SAr, ^oca, ^*Ti, 
*^Ci, ^^Fe, and ^^Ni (see also Nakamura et al. 2001 for 
details). The main purpose of the nucleosynthesis cal- 
culations is to see how much ^^Ni is produced at the 
explosion. For this purpose, inclusion of only a-nuclei is 
a good approximation because a-nuclei are the predom- 
inant yields of SNe. The equation of state takes into 
account gas, radiation. Coulomb interactions between 
ions and electrons, e~ — e"*" pair (Sugimoto & Nomoto 
1975), and phase transition (Nomoto 1982; Nomoto & 
Hashimoto 1988). To obtain many (i^kin — Mcj) rela- 
tions, we compute several hydrodynamical models with- 
out following nucleosynthesis and including only gas and 
radiation in the equation of state. The omitted physics 
in the equation of state, such as Coulomb interaction, 
mainly affects the result of nucleosynthesis and does not 
have much effect on i?kin and Mej. 

As the explosion mechanism of core-collapse SNe is not 
yet clear, we initiate the explosion as a thermal bomb 
(e.g., Nakamura et al. 2001). We put the thermal energy 
at Mr = 1.4 Mq in the nucleosynthesis calculations, as- 
suming that the 1.4 Mq neutron star is initially formed 
and the central remnant is treated as a point gravita- 
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Fig. 2. — Change in the velocity profile shows the shock prop- 
agation for the hydrodynamical model 13CO_2. The time since 
the explosion is shown. The mass cut between the ejecta and the 
fallback material is determined by whether the point exceeds the 
escape velocity or not. 
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Fig. 3. — Results of hydrodynamical calculations in the diagram 
of the kinetic energy of the ejecta (i?kin) a-nd the ejecta mass (Mej). 
Eltin is the final kinetic energy at infinity. For comparison, the lines 
oc e}I^ are also shown. 
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tional source. Here, is the mass coordinate from the 
center. There exist several ways to induce SN explosions, 
e.g., a kinetic piston (e.g., Woosley & Weaver 1995), but 
it is suggested that the results of nucleosynthesis are not 
sensitive to how energy is injected (Aufderheide et al. 
1991). However, note that Young & Fryer (2007) reports 
that the amount of fallback in low energy explosions de- 
pends on the method by which explosions arc indeuced. 
Generally, explosions by kinetic piston have less fallback 
because they tend to create stronger shocks and thus, 
the difference in the method also affects the yields of the 
nucleosynthesis. 

3.2. Bolometric Light Curve 

Bolometric LCs are calculated by using an LTE ra- 
diative transfer code (Iwamoto et al. 2000). This code 
assumes a gray atmosphere for the 7-ray transport. For 
the optical radiation transport, electron scattering and 
line opacities are taken into account. Electron number 
density is evaluated by solving the Saha equation. For 
simplicity, the line opacity is assumed to be a constant 
0.06 cm^ g~^. This value has been previously used for 
the explosion of CO stars (Maeda et al. 2003a). The gray 
7-ray opacity is set to be 0.027 cm^ g~^, which is known 
to be a good approximation (Axelrod 1980). Positrons 
emitted by the decay of ^^Co are assumed to be trapped 
in situ. To compare with the computed bolometric LCs, 
the observed bolometric LC of SN 2008ha is constructed 
as shown in the Appendix. 

4. RESULTS OF HYDRODYNAMICAL CALCULATIONS 
4.1. ii'kin ~ -^^cj Relation 

We calculate the hydrodynamics of the explosions and 
fallback for 25He, 40He, 13C0, 25CO, and 40CO with 
various input energies. In Figure [21 the change in the 
velocity profile shows the propagation of the shock wave. 
When the explosion energy i?kin is low, the inner part of 
the progenitor cannot overcome the gravitational poten- 
tial provided by the central remnant and thus falls back 
on to the remnant. In such cases, only the outer layers 
of the progenitor are ejected. The boundary between the 



fallback region and the ejecta is determined by whether 
the velocity of the region exceeds the escape velocity. 
We set a mass cut at this boundary to determine Mej. 
The expansion of the region above this mass cut eventu- 
ally becomes homologous. These homologous models are 
used for the calculations of the bolometric LCs. 

In Table [U we summarize i^kin and M^j for all our hy- 
drodynamical models and plot them in Figure [3] All the 
models with fallback are found to be on the line of either 

1 /3 

Afcj oc i?kin or Mej oc Ey.^^ . The reason why there exist 
two relations between ii^kin and M^j can be understood 
as due to a difference in the density structure. The den- 
sity structure of the progenitor affects the manner of the 
shock propagation, thus leading to the difference in the 
-E-kin -~ Mcj relation. 

Suppose the density structure of the progenitor is ex- 
pressed as p oc r~", where p is the density and r is the 
radius. Sedov (1959) showed that a shock wave is accel- 
erated when it is propagating along the density structure 
with a > 3 while it decelerates when propagating along 
the density structure with a < 3. This means that to 
achieve a certain velocity at a place with a density struc- 
ture of a < 3, more energy is required than the case of 
a > 3. As the escape velocity determines the boundary 
between the ejecta and the fallback region, the bound- 
ary is expected to be closer to the central remnant for 
the case for a > 3 than that of a < 3 if the same en- 
ergy is injected. Our hydrodynamical models show that 
if a > 3 at the boundary between the ejecta and the fall- 
back region, the results follow the relation Mcj oc Skin, 

1 /3 

and if a < 3, they follow A/oj oc E'j^jj^ . The exact phys- 
ical reason why Moj and i?kin are related as Mcj oc i?kin 

1 /3 

and Mcj oc E'j^jj^ is still unclear. We just treat it as an 
empirical relation in this paper and leave it as an open 
question. 

Figured] (left) shows the density structure of the pro- 
genitor model 25He. In Figure U (right), we plot hydro- 
dynamical models (Nos. 1-19) for 25He listed in Table 
[TJ The numbers attached to the points in the left panel 
of the figure are the same model numbers as in the right 
panel. The location of the model points indicates the 
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Fig. 4. — Density structure of 25He and the -Ekin ~ ^ej relation of 25He. Along the structure line, wc plot the points where the density 

slope is p oc (open squares). In the right panel, we plot hydrodynamical models (Nos. 1-19) for 25He listed in Table[T] The numbers 

attached to the points in the left panel are the same model numbers as in the right panel. The location of the model points indicates the 

mass cut of the hydrodynamical model. Looking at the density structure from outside, there is a small region where the density structure 

follows a < 3 which corresponds to model 3 with M^j ~ 0.10 Mq. Thus, the E^i^ — M^j relation of the models around this region follows 

1/3 1/3 
Mej OC Ejjijj . At M(,j ~ 0.33 — 0.81 Mq, the corresponding models 8-13 have a density structure with o < 3, thus following Af^j oc E^^^^ . 

At larger M^j, models 14-17 have a density structure with a > 3 and follow Af^j oc i?kin- The last two models, 18 and 19, have no fallback. 



TABLE 1 

Results of Hydrodynamical Calculations 



13CO 



25He 



25CO 



40He 



40CO 



No. 






Ekin 






A^ej 


Skin 




Ekin 




1 


0.00042 


0.035 


0.00035 


0.060 


0.0073 


0.056 


0.012 


0.10 


0.0038 


0.029 


2 


0.0012 


0.074 


0.00062 


0.091 


0.012 


0.093 


0.019 


0.14 


0.0050 


0.037 


3 


0.0015 


0.086 


0.0011 


0.12 


0.022 


0.17 


0.026 


0.19 


0.0072 


0.053 


4 


0.0018 


0.098 


0.0019 


0.16 


0.033 


0.24 


0.033 


0.23 


0.0094 


0.066 


5 


0.0022 


0.11 


0.0028 


0.21 


0.055 


0.39 


0.039 


0.29 


0.015 


0.10 


6 


0.0031 


0.14 


0.0038 


0.24 


0.074 


0.58 


0.054 


0.37 


0.027 


0.18 


7 


0.0042 


0.17 


0.0050 


0.29 


0.11 


0.78 


0.073 


0.51 


0.039 


0.25 


8 


0.0055 


0.22 


0.0078 


0.37 


0.14 


0.98 


0.11 


0.78 


0.066 


0.42 


9 


0.0071 


0.26 


0.011 


0.43 


0.19 


1.3 


0.15 


0.95 


0.093 


0.73 


10 


0.010 


0.32 


0.017 


0.51 


0.23 


1.5 


0.21 


1.2 


0.15 


1.3 


11 


0.013 


0.37 


0.022 


0.56 


0.30 


1.7 


0.35 


2.2 


0.21 


1.7 


12 


0.020 


0.45 


0.036 


0.69 


0.35 


1.9 


0.56 


4.1 


0.31 


2.5 


13 


0.028 


0.50 


0.045 


0.75 


0.43 


2.7 


0.86 


6.8 


0.39 


3.2 


14 


0.037 


0.54 


0.051 


0.83 


0.58 


2.8 


1.3 


9.7 


0.54 


4.4 


15 


0.047 


0.57 


0.060 


0.95 


0.74 


3.1 


1.7 


11 


0.66 


5.4 


16 


0.070 


0.64 


0.080 


1.6 


0.92 


3.2 






0.87 


7.1 


17 






0.12 


2.2 










1.0 


8.1 


18 






0.19 


4.4a 










1.3 


9.5 


19 






0.26 


4.4a. 










1.4 


10 


Note. — The units of E^iu and M^j are 


10^^ erg and Mq, respectively. 







No fallback 



mass cut of the hydrodynamical model. Looking at the 
density structure from outside, there is a small region 
where the density structure follows a < 3 which cor- 



portionality between i?kin and M^j was also shown by 
Nadezhin & Frank-Kamenetskii (1963) in the context of 
nova explosions. Thus, the ejecta velocity, which is scaled 



responds to model 3 with A/q ~ 0.10 Mq. Thus, the as w cx (^;ki„/Mej)^/^ would be a constant, being inde- 



^-kin ^ .^^oj relation of the models around this region fol 



1 /3 

lows Moj oc . 



At Mo 



0.33 - 0.81 Mn, the cor- 



responding models 8-13 have a density structure with 



1 /3 

a < 3, thus following Afej oc E^^l^ . At larger A/oj, 
els 14-17 have a density structure with a > 3 and follow 
Moj oc i^kin- The last two models, 18 and 19, have no 
fallback. 

For the extremely low i?kin, the mass cut in the explo- 
sion models lies in the outermost layer where the den- 
sity declines exponentially in all the progenitor models. 
Thus, the model sequences follow the relation i^kin oc M^j 



mod- 



pendent of ii^kin for each progenitor model. This means 
that however low Skin is, the ejecta velocity, i.e., the line 
velocities of the spectra, does not become lower than a 
certain asymptotic value. 



4.2. Rise Time 
Based on the i?ki 



Ejecta Velocity Relation 



Mej relation, we can construct a re- 



lation between observable quantities: the rise time (r) of 
the LC versus velocity (v) of the ejecta. The ejecta veloc- 
ity approximates the line velocities in the observed spec- 
tra. For each set of (i?kin, M^j), t is derived from the rela- 



at the low E],in hmit in the E],in - Moj plane. This pro- tion r oc k^/'^{M^:/ EmuY^^ (Arnett 1982), where k is the 
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Fig. 5. — Results of hydrodynamical calculations in the diagram 
of the rise time (t) and the ejecta velocity (v) defined by Equations 
((TJ and ((2]l. The attached numbers correspond to the numbers of 
the models shown in Table [l] The line on the right-hand side is 
the line Mgj = constant. 



total opacity and v is simply scaled as w oc (i?kin/A^oj)^^^- 
For simplicity, we assume that k, is constant for all the 
models. For illustration, we choose the proportional con- 
stants in T and v to match the typical values of Type la 
SNe, E^kin = 1-4 X 10^1 erg, Moj = 1.4 M©, r = 19.5 
days, and v — 9000 km s^^, and get the relations: 

jiM^y^\ (1) 

\^Skin/105iergy ^ ' ^ ^ 

. = 9000 f ^-;f;;> )^'\ms-^ (2) 

In Figure [5] we plot (t, v) for the model sequences in 
Table [T] and Figure [H As derived from Equations ([T]) 
and ([2]), the models with oc i?kin are located along 

1/2 1/2 

the line with r oc A/^j oc -Ekin ^^"-^ — constant, while 

1 /3 

the models with A/oj oc E^!^^ are located along the line 

1 /3 

with r ~ constant and v oc Mej oc E^^^ . For large -Bkin, 
all materials outside the proto-neutron star are ejected 
without fallback and Mcj = Mp^o — M^-cm is a constant, 
where Mp^o is the progenitor's pre-SN mass and Mrcm is 
the mass of the central remnant below the mass cut. As 
Moj is a constant, r and v follow the curve r^w oc Mcj — 
constant, as derived from Equations ([TJ and The 
black curve on the right side of Figure [5] shows the curve 
of Afoj = constant for model 25He. 

Our hydrodynamical models have a wide range of r 
as seen in Figure [S] and it implies that fallback SNe have 
a variety of LCs. However, v has only limited values for 
each progenitor and the progenitor of a fallback SN can 
be constrained not by its LC but by its photospheric 
velocity. In the next section, we constrain the progenitor 
of SN 2008ha mainly by using its photospheric velocity. 
As the observed line velocities around the maximum 
luminosity of SN 2008ha are 2000 km s-\ it is 
expected that 25CO, 40IIe, and 40CO would have too 
high photospheric velocities (w) to be consistent with 
those of SN 2008ha. 




velocity (cm 5^) 



Fig. 6. — Density structures of 13CO_2 at 2 days after the 
explosion. 



5. SN 2008HA 

We calculate bolometric LCs and the photospheric ve- 
locities for the models shown in Figure |3l We adopt 
the structure of the ejecta when it reaches the ho- 
mologous expansion in hydrodynamical calculations and 
assume that ^^Ni is uniformly mixed throughout the 
ejecta. Among all the models shown in Figure [31 we 
find that model 13CO_2 with (£;kin,Moj) = (1-2 x 
10^® erg, 0.074 Mq) is consistent with both the bolomet- 
ric LC and the photospheric velocity of SN 2008ha. Fall- 
back occurs in 13CO_2 and only the outermost layer of 
the progenitor is ejected. The density structure of the 
model is shown in Figure [6} Figure [71 shows that the cal- 
culated bolometric LC of 13CO_2 is in good agreement 
with the observed bolometric LC of SN 2008ha (see the 
Appendix). The '^^Ni mass ejected is assumed to be 0.003 
Mq- The rise time of 13C0_2 is 9.8 days. LCs from other 
explosion models of 13C0 are also shown for comparison 
in Figure [71 

In the nucleosynthesis calculation, the explosion of 
13C0.2 produces 0.15 M© of ^^Ni at Mr > 1.4 Mq. 
If we assume uniform mixing of ^^Ni at Mr > 1.4 Mq, 
the ejecta will contain 0.0086 Mq ^^Ni (Table ^. To 
reproduce the luminosity of SN 2008ha, 0.003 Mq of 
^^Ni needs to be contained in the ejecta. This implies 
that mixing due to the Rayleigh- Taylor instability (e.g., 
Hachisu ct al. 1991; Joggerst et al. 2009) or a jet (e.g., 
Maeda & Nomoto 2003b; Tominaga 2009) occurs in SN 
2008ha to bring ^^Ni to the outermost layer before the 
fallback. 

Figure [HI shows the photospheric velocities of 13CO_2 
compared with the observed line velocities of SN 2008ha 
(Figure 5 of F09). Among the line velocities shown in 
F09, we take Na i D and O i 7774 as good tracers of the 
photospheric velocity because these lines can be clearly 
distinguished and their line velocities are slower than 
other lines. The evolution of the photospheric velocity 
of 13CO_2 follows the line velocities of these tracers well. 
Thus, it is expected that our models will also be consis- 
tent with the observed spectra of SN 2008ha. Detailed 
synthetic spectra based on our model will be shown in a 
forthcoming paper (D. N. Sauer et al. 2010, in prepara- 
tion). 



6 



Moriya et al. 



TABLE 2 

Uniform Composition Based on the Result of the Nucleosynthesis of 13C0_2 



4He 


16o 


20Ne 


24Mg 


28 Si 


32s 


36 




44 Ti 


48 Cr 52 pg 


56Ni 


0.037 0.20 


0.32 


0.14 


0.066 


0.063 


0.034 


0.0065 


0.0045 


0.00051 


0.0017 0.011 


0.12 



Note. — The mass fraction of each element is shown. 




-11 L-U 1 1 1 1 ' ' 1 
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Fig. 7. — Bolometric LC of SN 2008ha and the calculated 
bolometric LCs for some explosion models of 13CO are shown. 
The bolometric LC of SN 2008ha is constructed as explained in 
the Appendix. Horizontal axis represents days since the explosion 
of the theoretical LCs of 13CO_2 and 13CO_5. The explosion date 
of 13CO_l is + 1 day in the figure. The kinetic energy and mass of 
ejecta of 13C0.2 are E^in = 1.2 X 10*8 erg and M^j = 0.074 Mq. 
The mass of ^^Ni is 0.003 Mq in 13C0.2. The rise time of 13C0.2 
is 9.8 days. E^in and M^j of 13CO_l and 13C0.5 are (Ekin, A^cj) = 
(4.2xlo4^erg,O.O35M0) and (2.2 X 10*8 erg, O.llM©), respectively. 

The model from 25CO whose photospheric velocity is 
shown in Figure [8] has £^kin — 3.3 x 10^^ erg and Mcj 
= 0.24 Mq (25C0_4). This model reproduce weh the 
bolometric LC of SN 2008ha. However, the photospheric 
velocity of this model is much higher than the line veloc- 
ities of SN 2008ha. This result is expected from Figure [5] 
because the ejecta velocities (v) of the explosion models 
of 25CO are too high for SN 2008ha (Figure [S|). On the 
other hand, the explosion model of 25He with the small- 
est energy we calculated (i^kin = 3.5 x 10'''' erg) still has 
a long rise time compared to SN 2008ha. Its LC is close 
to that of the explosion model 13CO_5 shown in Figure 
[3 which has almost the same r (Figure [5]). It is expected 
from Figures [3] and [5] that an explosion model of 25He 
with smaller energy could be consistent with SN 2008ha 
but the energy will have to be quite small (E'kin — 10^^ 
erg). Still, there remains a possibility that such an SN 
with a very small explosion energy could emerge as a 
result of the fallback. 

Although we assume a large amount of mass loss dur- 
ing the evolution of the progenitors, we do not assume 
the existence of a circumstellar matters due to a mass 
loss in our LC calculations. If the circumstellar matters 
are dense enough, SNe with fallback will be brightened by 
the interaction between the ejecta of SNe and the circum- 
stellar matters (Fryer et al. 2009). According to Fryer 
et al. (2009), this interaction could result in SNe having 
LCs with rising times and luminosities similar to those of 
SN 2008ha. Considering the fact that we do not see any 
evidence for the interaction in the spectra of SN 2008ha 




I ' ' ' ' 1 

10 20 30 40 50 

days since the explosion 



Fig. 8. — Photospheric velocities of 13CO_2 obtained from the 
bolometric LC calculations. These are compared with the observed 
line velocities of Na I D and O I 7774 from Figure 5 of F09, which are 
expected to be good tracers of the photosphere. Days are measured 
since the explosion date of each model. The rise time of 13CO_2 
is 9.8 days. We also show the explosion model of 25CO which has 
Ekin = 3.3 X 1049 erg and Afcj = 0.24 Mq. This model reproduces 
well the LC of SN 2008ha but the photospheric velocity is too high 
to be comparable to the line velocities of SN 2008ha. 

and there is no sudden drop in the tail of the LC of SN 
2008ha, which is expected in the interaction-powered SN 
models and is naturally explained by the nuclear decay 
of ^^Co, we think that SN 2008ha is likely to be pow- 
ered by the nuclear decay of ^^Ni and the circumstellar 
matter around the progenitor of SN 2008ha is so thin 
that the interaction does not become the major energy 
source of the LC. However, in conditions such as the mix- 
ing of the whole exploding star does not occur, almost 
all ^^Ni would fall back to the central remnant and only 
the interaction between the ejecta and the circumstellar 
matter would be the energy source to brighten the SN. 
We still do not have a clear observational SN with fall- 
back brightened by the interaction, but such SNe might 
be discovered in the future. 

6. DISCUSSION AND CONCLUSIONS 

We have found a fallback SN model (13CO_2) whose 
bolometric LC and photospheric velocity are both in 
good agreement with the observations of SN 2008ha. 
The ejecta of the model has very low explosion ener- 
gies and ejecta masses: (i?kin, M^}) = (1.2 x 10^^ erg, 
0.074 Mq). The explosion models from the progenitors 
of 25CO, 40CO and 40He are not in agreement with the 
observations of SN 2008ha because the photospheric ve- 
locities of these models are too high to be compatible 
with the observed line velocities of SN 2008ha (Figure 

. One might think that a model with sufficiently small 
-Ekin can have a low photospheric velocity. However, as 
shown in Figure [31 smaller i?kin leads to a larger fallback 
and thus smaller Mej as Moj oc i^kin- In this case, the ve- 
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Fig. 9. — Samples of the SED to construct the bolometric LC. The top left SED corresponds to the first point of the bolometric LC. 
Linear interpolations are made to fill the gaps in observations. The details of the interpolation are described in the Appendix. 



locity of the ejecta, v cx (E'kin/Afej) , will not become 
smaller. Alternatively, the explosion of 25He requires 
ii^kin — 10^^ erg to have a rise time similar to SN 2008ha. 

The difference between the successful and the unsuc- 
cessful progenitor models to explain SN 2008ha stems 
from the difference in the density structure. As discussed 
in Section 4.1, the density structure affects the propaga- 
tion of the shock wave and thus the relation between i^kin 
and Mcj . This means that whether a model can have the 
appropriate (i?kin,M3j) for SN 2008ha mainly depends 
on the density structure of the progenitor and not on the 
progenitor mass. Thus, the main-sequence mass of SN 
2008ha cannot be constrained only by hydrodynamical 
calculations and 13C0 would not be a unique progeni- 
tor candidate for SN 2008ha. Other progenitor models 
with appropriate density structures are also expected to 
reproduce the observational properties of SN 2008ha. 

The fact that faint SNe like SN 2008ha could emerge 
from the fallback of massive stars is also related to long 
gamma-ray bursts (LGRBs) without accompanied SNe 
(see also V09). LGRBs are thought to result from the 
death of massive stars and the fact that nearby LGRBs 
are accompanied by bright SNe (e.g., GRB 030329 and 
SN 2003dh, Hjorth et al. 2003; GRB 100316D and SN 
2010bh, Chornock et al. 2010) is one of the evidence for 
the scenario that relates LGRBs to the death of massive 
stars. However, some LGRBs are not accompanied by 
SNe even though they are close enough to be observed 
(e.g., GRB 060614 (Gehrels et al. 2006; Fynbo et al. 



2006; Delia Valle et al. 2006; Gal- Yam et al. 2006)) 
and this is one of the challenging problems of the the- 
ory of LGRBs. Several scenarios are proposed, e.g., a 
neutron star-white dwarf merger (e.g.. King et al. 2007) 
and a massive stellar death with a faint/dark SN (e.g., 
Tominaga et al. 2007). In this paper, we show that 
core-collapse SNe with fallback can be very faint and 
reproduce the observations of a faint SN (SN 2008ha). 
This supports the scenario that LGRBs without observed 
bright SNe are accompanied by very faint SNe with fall- 
back. In fact, theories of LGRBs like the coUapsar model 
(Woosley 1993) assume black hole formation with an ac- 
cretion disk and this picture is consistent with our fall- 
back SN model in the sense that most part of progenitors 
accretes to the central remnant and the central remnant 
becomes massive enough to be a black hole. This process 
could induce an LGRB, and it is possible that an LGRB 
was actually associated with SN 2008ha. We could have 
missed the LGRB because the jet of the LGRB was not 
directed to the Earth. In addition, our fallback SN model 
for SN 2008ha requires some mixing process to provide 
the ejecta with ^^Ni and the jet from the LGRB could 
have played a role in the mixing. 

Currently, several SNe are reported to have features 
similar to SN 2008ha. SN 2005E is one of the examples 
(Perets et al. 2010; see also Figure 7 of F09). As SN 
2005E is found far from the disk of the host galaxy, SN 
2005E is unlikely to be a core-collapse SN. However, its 
spectra do not show the characteristics of thermonuclear 
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Days since the explosion Days since the explosion 

Fig. 10. — Comparison of the bolometric LCs with those obtained in previous works. We apply our method to another SN 2002cx-Hke 
supernova SN 2005hk (PhiUips et al. 2007, Sahu et al. 2008) and a normal Type la supernova SN 2005cf (Wang et al. 2009). The rise 
times of SN 2005hk and SN 2005cf are assumed to be 15 days (Phillips et al. 2007) and 18 days (Wang et al. 2009), respectively. Our 
bolometric LCs tend to be a bit fainter around the maximum epoch but the decline rate and the luminosity at later epochs are in good 
agreement with other bolometric LCs obtained in previous works. 



explosions. Thus, SN 2005E is suggested to be a new 
type of stellar explosion (Perets et al. 2010). Perets et al. 
(2010) showed that the late phase spectra of SN 2005E 
contain strong emission lines of [Ca ii] AA7291, 7323 and 
suggested that those are the results of Ca enrichment in 
the ejecta. As this feature is similar to SN 2008ha (Figure 
7 of F09), they suggested that SN 2005E is related to 
SN 2008ha. However, the line velocities of SN 2008ha 
(^ 2000 km s~^) are considerably slower than SN 2005E 
(- 11,000 km s^i). Thus, it is not so obvious that SN 
2005E has the same origin as SN 2008ha. Kawabata et 
al. (2010) reported that Type lb SN 2005cz is another 
example of an SN that shows prominent Ca emissions and 
they related it to the core-collapse SN from a low mass (~ 
10 Mq) star. In this paper, we suggest that SN 2008ha 
is of core-collapse origin. Other SN 2002cx-like SNe, SN 
2002CX (Li et al. 2003; Jha et al. 2006) and SN 2005hk 
(Phillips et al. 2007; Sahu et al. 2008), are considered to 
be weak thermonuclear explosions while SN 2005E seems 
to be neither a core-collapse nor a typical thermonuclear 
explosion (Perets et al. 2010). Thus, SN 2002cx-like 
SNe might contain SNe with several kinds of origins and 
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Fig. 11.— Obtained bolometric LC of SN 2008ha. The rise time 
is set to 10 days in this plot. 



other criteria to classify them might be required. We 
need more samples of SN 2002cx-like SNe to clarify such 
criteria. See the supplementary information of Kawabata 
et al. (2010) for more intensive discussion for these Ca- 
rich SNe. 

We propose that SN 2008ha can be a core-collapse SN 
with fallback. However, it is not obvious that the fea- 
tures of the intermediate-mass elements that appeared 
in the earliest spectrum of SN 2008ha (Foley et al. 2010) 
can be synthesized by our model and this will be in- 
vestigated in a forthcoming paper (D. N. Sauer et al. 
2010, in preparation). Our model does not exclude the 
possibility that SN 2008ha is a thermonuclear explosion. 
However, simple estimates of Afoj by V09 and F09 are far 
below the Chandrasekhar mass limit (1.4 M©), which is 
required to ignite a thermonuclear explosion of a white 
dwarf. However, these estimates assume that the total 
opacity k is constant. It is possible that the effects of 
the opacity are large enough to create the appearance of 
a thermonuclear explosion similar to SN 2008ha. If this 
is the case, the opacity must become very low to match 
the fast rise of SN 2008ha as expected from the relation 
r oc K^''^(Mgj/£'kin)^/*. Other kinds of explosions such 
as Ma' SNe (e.g., Shen et al. 2010) or accretion-induced 
collapses (e.g., Darbha et al. 2010) might also be candi- 
dates but there still do not exist sufficient models that 
are consistent with SN 2008ha. 
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APPENDIX 

CONSTRUCTION OF BOLOMETRIC LIGHT CURVE 

We construct the bolometric LC of SN 2008ha by using the UBVRIJHK band LCs presented by F09. For each 
epoch, we derive the spectral energy distribution (SED) and obtain the bolometric luminosity by integrating the SED. 
To estimate the absolute magnitude of each band, we use the distance modulus /i = 31.64 mag and the color excess 
E{B — V) ^ 0.08 mag (F09). Reddening is corrected by following Cardelli et al. (1989). For each epoch, we use 
the observational data if they are available. If there is no observation of a band at a certain epoch, we deduce the 
magnitude of the band by linearly interpolating the observed magnitude of the nearest epochs. If no observational 
data are available before or after the observation as in the first few epochs when no RIK band data are available, we 
linearly interpolate data points on the SED. Finally, we integrate the SED by constructing trapezoids and triangles by 
connecting the UBVRIJHK data and the two edges. The two edges are chosen to be 1.25 x 10^^ Hz and 8.35 x 10^'' 
Hz, close to the frequencies of the K band and U band. Some samples of the SED are shown in Figure [9l We calculate 
the bolometric LCs of SN 2005hk (PhiUips et al. 2007, Sahu et al. 2008) and SN 2005cf (Wang et al. 2009) to check our 
method (Figure [TU]) . Our bolometric LCs are in a good agreement with the bolometric LCs obtained in the previous 
studies. 

The bolometric LC of SN 2008ha obtained by our method is shown in Figure [TTJ As the first point is constructed 
by only the JH data, we assume that the actual luminosity at this epoch is higher. The rise time of our bolometric 
LC is consistent with that of the quasi-bolometric LC constructed by F09. 
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